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ABSTRACT 

In this letter we carry out the first systematic investigation of the expected gravitational wave 
(GW) background generated by supermassive black hole (SMBH) binaries in the nHz fre- 
quency band accessible to pulsar timing arrays (PTAs). We take from the literature several 
estimates of the redshift dependent galaxy mass function and of the fraction of close galaxy 
pairs to derive a wide range of galaxy merger rates. We then exploit empirical black hole- 
host relations to populate merging galaxies with SMBHs. The result of our procedure is a 
collection of a large number of phenomenological SMBH binary merger rates consistent with 
current observational constraints on the galaxy assembly at z < 1.5. For each merger rate we 
compute the associated GW signal, eventually producing a large set of estimates of the nHz 
GW background that we use to infer confidence intervals of its expected amplitude. When 
considering the most recent SMBH-host relations, accounting for ultra-massive black holes 
in brightest cluster galaxies, we find that the nominal la interval of the expected GW signal 
is only a factor of 3-to-10 below current PTA limits, implying a non negligible chance of 
detection in the next few years. 

Key words: black hole physics - gravitational waves - pulsars: general - galaxies: evolution 



1 INTRODUCTION 

Precision timing of an array of millisecond pulsars (PTA) pro- 
vides a unique opportunity to get the very first low-frequency 
gravitational wave (G W) detection. The Eu ropean Pulsar Tim- 
ing Array (EPTA, iFerdman et all |2010|). the Parkes Pul- 



sar Timing Array (PPTA, Manchester et alj 1201 2|) and the 



North Americ an Nanohertz Ob servatory for Gravitational Waves 
(NANOGrav, Ijenet et alj|2009h , joi ning together in th e Interna- 
tional Pulsar Timing Array (IPTA, iHobbs et alj|201(]|) . are con- 
stantly improving their sensitivity in the frequency range of ~ 
10 -9 — 1CP 6 Hz. Inspiralling supermassive black hole (SMBH) 
binaries populating merging galaxies throughout the Universe 
are expected t o generate the dominant signal in this frequency 
band (see, e.g. iRaiagopal & Romaij|l995t IJaffe & Backerlbdol ; 
IWvithe & Loebl l2003; Sesana et al. 2008]). Generally speaking, the 
expected amplitude of the signal depends on the pace at which 
SMBH binary mergers occur along cosmic history, and on their 
typical masses. Both quantities are poorly determined observa- 
tionally, allowing for a wide range of GW signal amplitudes. 
Theoretical models of SMBH evolution within the standard hi- 
erarchical framework of galaxy formation indicat e a typical GW 
strain amplitude A ~ 10~ 15 at / = 1/yr dWvithe & Loebl 
l2003l : ISesana et al.ll2008l:lRavi et alj|2012l) , with an uncertainty of 
«0.5dex i Sesana et al. 20081) . However, a recent investigation by 
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iMc Williams et al.1 d2012l) . based on a phenomenological model in 
which the low redshift massive galaxy assembly is driven by merg- 
ers only, predicts a higher background with a fiducial amplitude 
A ~ 6 x 10~ 15 . Though useful, all the aforementioned models em- 
ploy specific recipes for the galaxy assembly and/or for the growth 
of SMBHs, and a systematic investigation of the possible range of 
signals compatible with observational uncertainties is still missing. 
This is a particularly important issue to assess at this point for two 
reasons: (i) the best limit placed by PTAs on the GW back ground 
amplitude is A — 6 x 10 -15 dvan Haasteren et all 1201 ll) . close 
to theoretical predictions; (ii) SMBHs in brightest cluster galax- 
ies (BC Gs) were recently discovered to be more massive than ex- 
pected ( Hlavace k-Larrondo et alj|2012l), resulting in a revisio n to 
the established SMBH-host relations dMcConnell & Mall2012h that 
might push the expected GW background level closer to current 
upper limits, implying possible detection in the next few years. 

As part of the common effo rt of the EPTA collaborati on to 
detect GWs with pulsar timing dvan Haasteren et al. I I20TT1) . we 
present here the first systematic investigation of the range of 
GW signal amplitudes consistent with observationally based es- 
timates of the SMBH assembly in the low redshift Universe. The 
manuscript is organized as follows. In Section 2 we describe our 
model for generating the GW background and we test it against 
a range of observational constrains. We present and discuss our 
main results in Section 3, and draw our conclusions in Section 4. 
Throughout the paper we assume a concordance A-CDM universe 
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with Qm = 0.27, Q\ — 0.73 and h — 0.7. Unless otherwise 
specified, we use geometric units where G — c — 1. 



2 BUILDING THE GW BACKGROUND FROM 
ASTROPHYSICAL OBSERVABLES 

2.1 Mathematical description of the GW background 

Consider a cosmological population of merging SMBH binaries. 
Each merging pair is characterized by the masses of the two holes 
Mm i > M m , 2, defining the mass ratio q, = M. ,2/M.,fl Follow- 
ing [Sesanaetal ] j2008l) (see also |Phinnevll200ll) . the characteristic 
amplitude h c of the GW signal generated by such population is 
given by 

j3 ~ 1 dE gv (M) 



K(f) 



■Kf 2 



dzdM,,idq, 



d 3 n 



dzdMm idq,, 1 + z din f r 
(1) 

Here, the ene rgy emitted per logarithmic frequency interval is 
dThomell 19871) 



dEgw _ k 2/3 s/3 2/3 
dln/ r 3 /r : 



(2) 



where we assumed circular binaries driven by GW emission onl)[f], 
M = (M.,iM., 2 ) 3/5 /(-M»,i + M., 2 ) 1/5 is the chirp mass of 
the binary and f r — (1 + z)f = is the GW rest frame fre- 
quency, which is twice the binary Keplerian frequency. The quan- 
tity d 3 nj {dzdMm , i dq, ) represents the differential merger rate den- 
sity (i.e. number of mergers per comoving volume) of SMBH bi- 
naries per unit redshift, mass and mass ratio. For convenience, we 
decided to keep it a function of M m ,i and q, inst ead of M only. 

It is straightforward to show dPhinne y|200ll) that the predicted 
characteristic amplitude scales as oc ,f~ 2 ^ 3 , with a normalisation 
that depends on the details of th e merging binary p opulation, and is 
usually represented as (see, e.g. Jjenet et al.ll2006l) : 



ho(f) = A 



f 



yr" 



-2/3 



(3) 



where A is a model dependent constant that represents the ampli- 
tude of the signal at the reference frequency / = lyr~ . Since ob- 
servat ional limits on the GW background are usual l y given in terms 
of A l ljenet et alj|2006l : Ivan Haasteren et alJ|2oTH : IPemorest et al.l 
I2012T) . we keep the parametrization given by equation {3]l in this 
letter, and we investigate the range of A predicted by phenomeno- 
logical models of the SMBH assembly based on observations. 

Before proceeding, it is worth mentioning that the shortcom- 
ings of equation (TJ in catching the relevant features of the GW 
signal emitted by a realistic population of quasi monochromatic 
sources were ext ensively investigated by dSesana et alj2008l . l2009l ; 
iRavi et alj|2012h . Although equation Q} fails in describing small 
number statistics effects and the intrinsic non Gaussianity of the 
signal, it is sufficient to describe its expected overall amplitude, 
which is our main interest here. In a companion paper (Sesana in 
preparation) we will carry out a more systematic study of all the 
relevant signal features, including statistics of expected individu- 
ally resolvable sources. 

1 to avoid confusion, we denote the masses and mass ratio of the SMBH 
binary as Mm.i, Mm,2,Qm, whereas plain M and q are used for galaxies 

2 We keep this assumption throughout the paper. Eccentricity, together 
with other physical mechanisms (stellar scattering, gas torques) driving the 
binaries can modify the form of dE gvi /d In f r given by equation (2)- We 
defer the investigation of these issues to future work. 



2.2 Determination of the SMBH binary merger rate 

Since the energy emitted per logarithmic frequency interval is fixed 
by General Relativity (in the approximation of circular GW driven 
binaries), the typical background strength A depends on the SMBH 
binary differential merge r rate only. In contrast to our past work 
dSesana et alj|20 08, 2009). we take here an observational approach 
to determine d 3 n/dzdMm,idq,. We proceed in two steps: (i) we de- 
termine from observations the galaxy merger rate d 3 riG jdzdMdq 
(in a merging galaxy pair, M and q < 1 are the mass of the pri- 
mary galaxy and the mass ratio respectively), and (ii) we populate 
merging galaxies with SMBHs according to empirical black hole 
mass-galaxy host relations found in the literature. 



2.2.1 Galaxy merger rate 

The galaxy differential merger rate can be written as 

d 3 n G _ </>(M, z) T(z, M, q) dt r 
dzdMdq ~ Mm 10 r(z,M,q) ~dz' 



(4) 



Here, 4>(M, z) = (dn/dlogM) z is the galaxy mass function mea- 
sured at redshift z; T(M,q,z) — [df /dq)M,z is the differential 
fraction of galaxies with mass M at redshift z paired with a sec- 
ondary galaxy having a mass ratio in the range q, q+Sq; t(z, M, q) 
is the typical merger timescale for a galaxy pair with a given M 
and q at a given z\ and dt r /dz converts a proper time rate into a 
redshift rate and is given by standard cosmology. The reason for 
writing equation © is that (f> and T can be directly measured from 
observations, whereas r can be inferred by detailed numerical sim- 
ulations of galaxy mergers, as discussed below. 

We t ake three different galaxy stellar mass functions from the 
literature dBorch et alj|2006l ; iDrorv et alj|2009t lllbert et alj|2010h 
and match them with the local mass function dBell et alj|2003l) . to 
obtain three fiducial <f> z (M). To each fiducial mass function we 
add an upper and a lower limit accounting for the errors given by 
the authors on the function best fit parameters, plus an additional 
0.1 dex systematic error due to uncertainties in the determination 
of the galaxy masses, for a total of 9 galaxy mass functions. For 
all mass functions we separate early type and late type galaxies. 
We restrict our calculation to z < 1.3 and M > 1O 1O M0, since 
these are the systems contributing the largest fraction of the GW 
signal. By extrapolating our calculations we found that merging 
pairs residing in galaxies with M < 10 1 Mq or at z > 1.3 can 
contribute at most 5% to the signal amplitude, and can be safely 
neglected. 

We consider fo ur studies of the evolution of the 
galaxy pair fraction dBundv et al.l [20091: Ide Ravel et all 120091: 
iLopez-Saniuan etld] 120121 : IXu et al.l |2012|) . Pair fractions are 
usually integrated over some range of q and are given in different 
mass bins in the form ,f(z) = /o(l + z) 7 . A good proxy for the 
observed pair q distribution is df /dq cc q^ 1 . Given f(z), we can 
therefore simply write df/dq(z) — —f(z)/(qlnq m ), where q m 
is the minimum mass ratio selected in counting pairs. Each author 
applies different criteria as for q m , mass and redshift range, and the 
maximum projected distance d ma x below which two galaxies are 
considered a bound pair, as detailed in table [TJ Also in this case, 
for each of the 4 fiducial models we consider an upper and a lower 
limit taking into account the the errors in the best fit parameters 
fo, 7, as reported by the authors, to get a total of 12 pair fraction 
models. When necessary, we extrapolate the pair fraction estimates 
to cover the full mass and redshift range of interest (z < 1.3 and 
M > 10 10 Mo). When pair counting for different galaxy types are 
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Paper 


</m 


M min [M ] 


<2 ma x[kpc] 


gal. type 


Bundy et al. 2009 


0.25 


10 10 


20 


yes 


de Ravel et al. 2009 


0.25 


10 9.5 


100 


no 


Lopez et al. 2012 


0.25 


10 11 


30 


yes 


Xu et al. 2012 


0.4 


lfJ 9.4 


20 


no 



Table 1. Overview of the pair fraction selection performed in the paper used 
in this work. See text for details. 



available, we apply them to the corresponding galaxy type mass 
function, otherwise, w e assume the same p a ir frac tion for early 
and late type galaxies. lLopez-Saniuan et"al I J2012I) also provide 
pair fractions for 'minor mergers', i.e., for 0.25 > q > 0.1. 
We checked that including those in our calculation enhances the 
background by a factor 0.06dex (< 15%) at most. 

Galaxy merger times cales r were carefully estimated by 
IKitzbichler & White! |20o3) us ing mock catalogues of galaxy pairs 
in the Millennium simulation JSpringel et aT]|2005h . In their equa- 
tion (10) they provide the average merger timescale as a function 
of M, z, and projected distance d p . We complemented their equa- 
tion (10) with a ~ g -0 ' 3 dependence extracted by fitting the re- 
sults of a set of full hyd r odyna mical simulation of galaxy mergers 
presented bv lLotzetall feoich. In doing this, we noticed that the 
merger timescales given by Lotz et al 



l20ld) are a factor of two 



shorter than those given by IKitzbichler & White! J2OO8I) ; we there- 
fore adopted two different normalizations to get a 'fast' and a 'slow' 
merger scenario. 

We interpolate all the measured <f>, T , r on a fine 3-D grid in 
(z,M,q), to numerically obtain 9 x 12 x 2 = 216 differential 
galaxy merger rates. Note that typical values of r are of the order 
of a Gyrs, therefore, the merger rate at a given (z, M, q) point in 
the grid is obtained by evaluating 4> and T at (z + 5z, M, q), where 
Sz is the redshift delay corresponding to the merging time r. Note 
that by doing this, we implicitly assume that all SMBH binaries 
coalesce instantaneously at the merger time of their hosts. 

2.2.2 Black hole-host relations 

We assign to each merging galaxy pair SMBHs with masses drawn 
from 9 different SMBH-galaxy relations found in the literature (see 
tableO. We write them in the form 



log 1() A/. = a + @\og w X, 



(5) 



where X = {cr/200km s"\ Lt/W^L^Q or A'/ b uigc/10 11 M o }, 
being a the stellar velocity dispersion of the galaxy bulge, L; its 
mid-infrared luminosity, and Mb u i gc its stellar mass. Each relation 
is characterized by an intrinsic scatter e. a, (3,e are listed in table[2] 
The relations link Al, to the bulge properties, whereas our galaxy 
merger rates are function of the total stellar mass. We derive the 
bulge mass of each galaxy by multiplying the total stellar mass by 
a factor /bulge- We assume /bulge — 1 for 3-11 early type galaxies 
with M > 10 n M o , declining to f bulgc = 0.5 at M = 10 10 M Q , 
whereas we assign a random /bulge in the range 0.1-0.3 to late type 
(i.e., disk dominated) galaxies. Although this prescription is some- 
what arbitrary, we found that the results are almost independent of 
/bulge as long as massive early type galaxies retain a bulge frac- 
tion of order unity, which is a well established observational fact. 
Given the bulge mass, we est i mate L j by inverting the A'ibuigc — Li 
relation given by ISani et al. and we compute a by fitting 

a broken power-law to the z = a — A'/buigc data presented by 



Ptipcr 


v 


(V 


a 
P 


c 


Harine & Rix (2004) 


A/bulge 


8.2 


1.12 


0.30 


Sanietal. (2011) 


A^bulgc 


8.2 


0.79 


0.37 


Beifiori et al. (2012) 


-^bulge 


7.84 


0.91 


0.46 


McConnell & Ma (2012) 


■^bulge 


8.46 


1.05 


0.34 


Sanietal. (2011) 


Li 


8.19 


0.93 


0.38 


Gultekin et al. (2009) 


(7 


8.23 


3.96 


0.31 


Graham et al. (2011) 


(J 


8.13 


5.13 


0.32 


Beifiori etal. (2012) 


(7 


7.99 


4.42 


0.33 


McConnell & Ma (2012) 


a 


8.33 


5.57 


0.40 



Table 2. List of parameters a, fj and e. See text for details. 
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Figure 1. Upper plot: local SMBH mass function. Thick and thin solid 
black lines enclose the areas corresponding to 68% and 95% confidence 
levels given b y our models. Colored s haded areas are SMBH mass functions 
estima t ed by iMarconi et all {2004); Shan kar et al] 12004}) ; iHopkins et alj 
<2007l) ; iTundo et alj fe007l) . Lower plot: redshift evolution of the total 
SMBH mass density. Thick and thin solid black lines have the same mean - 
ing as in the upper plot. Red and green dots are from lZhang et all fe012l) . 
blue dots are from IHopkins et all j2007l) . cyan and yellow lines are from 
iMerloni & Heinz! |2008|) , magenta lines bracket the la uncertainty given 
bv lShankar et alj 120041) . a nd the thick ma genta line is the estimated uncer- 
tainty range at z = from Shankar (2009). 

iRobertson et alj d2006t) . When converting Mb u i ge into a we apply 
a multiplication factor (1 + z) a ' :i , to account for the observational 
fact that galaxies of a given mass at higher redshift are more con- 
centrated and have larger velocity dis persions than galaxies of th e 
same mass at lower redshift (see, e.g. lLopez-Saniuan et al.ll2012l ). 
Having derived A/b u i ge , Li and a, we can populate galaxies with 
SMBHs. We apply a further (1 + z) 0,3 correction to the scaling 
relations involving Afbuigc and Li. This redshift dependence im- 
proves the match of the SMBH mass density (pbh) redshift evolu- 
tion given by our models with other estimates found the literature 
(see figure O, and we checked that our results are basically inde- 
pendent of it. 

We assign to each merger remnant a total bulge mass equal 
to the sum of the total stellar masses of the merging systems, i.e., 
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/bulge = 1. Furthermore, we correlate the masses of the merging 
SMBHs either to the properties of the two merging galaxies or to 
those of the merger remnant, fol lowing the scheme described in 
Section 2.2 of lSesana et all j2009h . This gives us three slightly dif- 
ferent mass estimations for the SMBHs forming the binary for each 
adopted scaling relation. 

We combine the 9x3 = 27 different ways to populate 
the merging galaxies with SMBHs together with the 216 galaxy 
merger rates to obtain 5832 different SMBH binary merger rates 
d 3 n/dzdM, : idq m , consistent with current observations of the evo- 
lution of the galaxy mass function and pair fractions at z < 1.3 
and M > 1O 1O M and with the empirical SMBH-host relations 
published in the literature. We give equal credit to each model, 
and we generate 5832 GW signals, sufficient to place reasonable 
confidence levels for the expected amplitude according to current 
observational constraints. Our approach is modular in nature, and 
it is straightforward to expand the range of model to include new 
estimates of all the quantities involved. 



2.3 Validation of the models 

Although the evolution of the SMBH masses is not followed self- 
consistently in our models, in figure [TJ we validate them by com- 
paring the local SMBH mass function and the redshift evolution 
of the total SMBH density with several estimates found in the lit- 
erature. We also checked that the predicted range of galaxy and 
SMBH merger rates as a function of mass and redshift are broadly 
consistent (though with a large scatter) with those derived from 
our previous models con structed on top of the Millennium Simula- 
tionj Sesana et alj|2009l) or exploiting semianalytical merger trees 
dSesana et al .Il2008h . In the latter approach we evolve the SMBH 
population self-consistently. In figure [TJ we show the nominal la 
and 2a confidence levels (i.e. the range in which 68% and 95% 
of our models are contained) of the estimated local SMBH mass 
function and mass density as a function of z. The agreement with 
independent results published in the literature is excellent. We no- 
tice that we allow for slightly larger values of b oth quantities with 
respec t to published results. This is because the lMcConnell & Mai 
d2012l) scaling relations, that include the recently measured ultra- 
massive SMBHs in BCGs, predict SMBH masses which are 0.2-to- 
0.4dex larger than previous estimates at the high mass end. Those 
models will result in larger amplitude of the GW signal, which 
might be soon directly tested with PTA observations. 



3 RESULTS 

Our main result is shown in figure [2] where we plot confidence 
levels on the GW characteristic amplitude given by our models. 
When considering the whole set of models (upper left panel), the 
68% confidence region lies in the range 3.3 x 10~ 16 < A < 
1.3 x 1CP 15 , corresponding to a factor of 4 uncertainty in the 
GW signal. The 99.7% region extends much further, in the range 
1.1 x KT 16 < A < 4.2 x 1CT 15 , corresponding to a factor 
« 40 uncertainty. Note that this l atter upper bound is only a factor 
1.5 below the best limit placed by van Haa steren et alj J201 ll) . Our 
'democratic' approach to the problem gives the same weight to all 
the models. One can argue that models featuring the best estimates 
of the galaxy mass function and pair counts, should be considered 
more robust than those constructed using the upper or lower limits 
for the same quantities (see Section 2.2.1). If we restrict to 'fiducial 
models only', the scatter is mildly reduced, and the 68% and 99.7% 
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Figure 2. Characteristic amplitude of the GW signal. Shaded areas repre- 
sent the 68%, 95% and 99.7% (nominally la, 2a, 3a) confidence levels 
given by our models. In each panel , the b lack asterisk marks the best cur- 
rent limit from I van Haasteren et alj 4201 lh . Sha ded areas in the upper lef t 
panel refer to the 95% confidence level given bvlMc Williams et alj <20l3) 
(red) and the uncertainty range estimated bv lSesana*et*alTl2008h . See text 
for discussion. 




14 



Figure 3. Normalized distributions of the expected GW amplitude A at / = 
lyr -1 . Black solid line, all models; green dot-dashed line, fiducial models 
only; red short-dashed line, models antecedent SMBH measurements in 
BCGs; blue long-dashed, models including SMBH measurements in BCGs. 
The shaded area marks the region excluded by current PTA limits, whereas 
the solid dotted line represent what can be achieved by timing 20 pulsars at 
100ns rms precision for 10 years. 



confidence levels are set in the range 3.8 x 10 16 < A < 1.1 x 
10~ 15 and 1.7 x 10~ 16 < A < 2.2 x 10~ 15 respectively (upper 
right panel). Things become much more interesting if we consider 
only the SMBH-host relations updated to inclu de the recent mea- 
surements of ultra-massive black hole in BCGs dMcConnell & Mai 
|2012|) . As expected, the signal is boosted-up, bringing the 68% and 
99.7% confidence intervals to 5.6 x 10~ 16 < A < 2.0 x 10~ 15 
and 2.4 x 10~ 16 < A < 5.7 x 10~ 15 respectively (lower right 
panel), a factor ~ 2 larger then models featuring previous estimates 
of the SMBH-host relations (lower left panel). Although obtained 
with a completely different procedure, our confidence intervals 
are generally consi stent with the estimated sig nal range given by 
dSesana et alj|2008l) . whereas recent results by iMc Williams et"alj 
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d2012l) are marginally consistent (at a 3a level) with our findings. 
This is not surprising, since their purely merger driven SMBH evo- 
lution naturally produces the highest possible signal for a given 
SMBH mass function. Within a year, IPTA observations will there- 
fore be able to test the whole amplitude range predicted by this 
scenario (van Haasteren, private communication) . In figure [5] we 
plot the normalized distributions of A given by all our models. The 
overall distribution (solid-black line) has a neat gaussian shape, 
in agreement with the central limit theorem. This is the sign that 
none of the ingredients of our model (galaxy mass function and 
pair fraction, coalescence time, SMBH-host relation) plays a par- 
ticularly dominant role in determining the signal amplitude. The 
shaded area, marking the region of A excluded by current limits, 
already overlaps with the long tails of the distributions. Roughly 
speaking, the maximum A detectable by a PT A with a signal-to- 
noise ratio of 5 is given by l lSesana et alj|2008h 



A « 8 x 10" 



6t 



ioons viooy \2o) v.5yr 



-1/2 



-5/3 



(6) 

where 8t ms is the rms residual of each individual measurement 
(assumed to be the same for each pulsar), N r is the number of 
measured residuals for each pulsar, N p is the number of pulsars in 
the array and T b s is the duration of the experiment. Observations 
of 20 pulsars at 100ns rms precision for 10 years will allow to detect 
a signal of A w 2 x 10 -16 , which encompasses more than 95% of 
the models presented here (dotted vertical line in figure[3]l. 



4 CONCLUSIONS 

We presented the first systematic investigation of the GW back- 
ground generated by a cosmological population of SMBH binaries 
in the nHz frequency band, relevant to PTAs. We generated a grand 
total of 5832 SMBH binary merger rates, by coupling several ob- 
served galaxy mass functions and pair counts to phenomenological 
SMBH-host relations, and assuming merger timescale prescriptions 
derived by detailed hydrodynamical simulations of galaxy merg- 
ers. By construction, our models are consistent with observational 
constraints, and produce SMBH mass functions and mass density 
evolutions with redshift consistent with several independent esti- 
mates found in the literature. When considering all models, we find 
characteristic amplitudes of the signal at a frequency lyr -1 in the 
range 3.3 x 10~ 16 < A < 1.3 X 10~ 15 at 68% confidence. How- 
ever, recent measurement of ultra-massive black holes in BCGs led 
to a revision of the SMBH-host relations, rising their high-mass 
end by 0.2-0.4dex. Models based on these latter results predict 

and 

at 68% and 99.7% confidence 



amplitudes in the range 5.6 x 10 < A < 2.0 x 10 



2.4 x 10" 



< A < 5.7 x 10" 



respe ctively. Our results are broadly consistent with our previous 
work dSesana et al.ll2008|) and marginal l y cons istent (at a 3a level) 
with recent work by iMc Williams et al.l HoT3> . We predict a nom- 
inal 3cr upper limit to the signal close to current limits placed by 
PTAs. If our models are correct, with an improvement of a factor 
of only three on current limits, there is a non negligible chance to 
make the first ever direct GW detection. Even a negative result will 
nevertheless allow us to constrain the assembly of the most mas- 
sive galaxies at low redshift and how do they correlate with their 
hosts, turning PTAs into useful astrophysical probes. Looking fur- 
ther ahead, the timing of 20 pulsars at 100ns rms precision for 10 
years (considered a feasible long term goal with current PTAs) is 
almost guaranteed to detect the GW background. The first direct 
GW detection, might not be so far off in the future. 
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